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Abstract: In nonequilibrium thermodynamics macroscopic entropy creation plays an im- 
portant role. Here we study, from various viewpoints, its relation with the phase space 
contraction, which has been recently proposed as an apparently alternative quantity. 
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1. Thermostats, entropy creation and phase space contraction 

Studying stationary states of mechanical systems in interaction with thermostats the 
latter are often modeled by systems of particles subject to anholonomic constraints. The 
equations of motion take the form x = fE_{x) where x is a point in phase space and E_ are 
parameters controlling the size of the acting nonconservative forces. 

def 

It has appeared natural to define entropy creation rate as the divergence o(x) = — 
J2j d Xj fE,j{x). There are other natural definitions of entropy creation rate and here we 
study their relation with the above phase space divergence. The aim is to find such a 
relationship for a system that can be considered to be described by a continuum evolving 
under macroscopic equations. 




Fig.l Reservoirs occupy finite regions outside Cq, e.g. sectors C; C i — 1,2 . . .. Their particles are constrained to have a 
total kinetic energy K% constant, by suitable forces, so that the reservoirs "temperatures" Ti are well defined. 

To be concrete and in rather ample generality we imagine a system C of particles en- 
closed in a container, also called Co, with elastic boundary conditions surrounded by a few 
thermostats which consist of particles, all of unit mass for simplicity, interacting with the 
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system via short range interactions, through a portion diCo of the surface of Co, and subject 
to the constraint that the total kinetic energy of the iVj particles in the i-th thermostat is 
Ki = \X_ i = ^Nik B Ti. A symbolic illustration is in Fig.l. The equations of motion will be 

lo = -% (t/o(2C ) + E w oA2L ,K j )) + E(X ), 

X t = -dx t (U t (2Li) + Wo.&Xj)) - ajQ (1) 

with «j such that Ki is a constant. Here W 0ti is the interaction potential between particles in 
Ci and in Co, while Uo, Ui are the internal energies of the particles in Co,C« respectively. We 
imagine that the energies Woj,Uj are due to smooth translation invariant pair potentials; 
repulsion from the boundaries of the containers will be elastic reflection. It is assumed, 
in Eq.(l) that there is no direct interaction between different thermostats: their particles 
interact directly only with the ones in C . Here E_{Ko) denotes possibly present external 
positional forces stirring the particles in C . 

Since the work per unit time that particles outside the thermostat (hence in C ) exercise 
on the particles in it, is Qi d = — dy_ Wm (X . X { ) • X ; , and it can be interpreted as the "amount 
of heat Qi entering" the thermostat Ci, energy conservation yields 



j t {\x* + Ui) = U, = -a t X- + (2) 

and the contraints on the thermostats kinetic energies give ctj = 3 ^'~^ . 

Set x = (2Li, 2Li)i=o,.. an d write Eq.(l) as x = /b(x). The divergence a(x) = 
~ J2j d Xj fE,j(x) of the equations of motion in phase space is readily computed from Eq.(l), 
see also [1], and is a(x) = X^>o7pr + R or 



a(x) = e{x) + R, with 

e(*) = Eor (3) 

where R = J2i>o]^T~ : ^ ere anc ^ fcfr"'fci^" snom d really be multiplied by (1 — 3^-): this 
neglection of Oi^Nf 1 ) is made just to simplify the formulae as, in any event, we shall be 
interested in cases in which N , Ni 3> 1. 
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Remark: (i) e(x) = Z)i>Ofcfr" can ^ e cane d naturally the entropy creation rate and, therefore, 
Eq. (3) has a physical meaning: entropy is created at the boundary of the system. Creation 
really takes place where the walls get in contact with the thermostats, where the tempera- 
tures Ti are defined. 

(ii) Note that if particles in C were also subject to an isokinetic constraint |(2l ) 2 = 
^N k B T phase space contraction would simply be changed by the addition of ^jl. with Q 
being the work done per unit time by the thermostats in d, i > 0, on particles in C ; also R 
will contain an extra term proportional to Uq. 

(iii) The divergence a(x) is different from the entropy creation rate e(x). Their difference is 
a "total time derivative", see Eq.(3), therefore there is a relation between the time averages 
a d = \ f\ cr(Stx) dt and a d = \ fix. s(S t x) dt, namely 



o = a + ^(R(Sr X ) - R(S-rx)). (4) 

This means that the observables a and ao will have the same distribution with respect to 
any stationary distribution in the limit r — > oo if R is a bounded function (as in our case). 
More general and "singular" interaction potentials could be considered and would lead to 
essentially equivalent conclusions, following the results in [2]. 

(iv) Note that also phase space contraction of a system in contact with isokinetic thermostats 
has a precise physical meaning as it equals minus the sum of the dimensionless free energy 
creation rates —-r^r + 7%- of the thermostats. 

2. Macroscopic fluids 

The above analysis shows that the two notions of entropy creation rate e(x) in Eq.(3) 
and of phase space contraction a(x) are related but different. 

They have in particular the same stationary average, as they differ by a total derivative R. 
As a consequence, not only the averages of a and e are equal but also that the fluctuations 
of their finite time averages, i.e. of a and ao in Eq.(4), are the same in the limit r — > oo: 
so that properties known for the fluctuations of a imply corresponding properties for the 
fluctuations of the physically meaningful entropy creation rate e. This is relevant because, 
in the literature, several results have been derived concerning the fluctuations of the time 
averages of the phase space contraction, see [4]. 
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Therefore it is of some interest to see what the above mechanical notion of entropy 
creation rate becomes in a system which can be considered as a continuum in a stationary 
state and in local equilibrium. In fact for such a system an independent definition of entropy 
creation is classical, [3]. We check that the two notions coincide up to a total derivative. 

Consider, in C , a system of particles which can be regarded as a Navier-Stokes (NS) 
continuum in a stationary state and in contact with fixed walls on which, at each boundary 
point | G dCo, temperature is prescribed at a value T(£) because the surface element ds$ is 
in contact with a thermostat (as idealized in Fig.l) whose internal potential energy per unit 
surface is U ext (£). Then the entropy creation rate according to Eq.(3) will be 



where Q(£) is the amount of work per unit time and unit surface that the fluid performs on 
the thermostat in contact with the surface element ds$, while phase space contraction a will 



There are no complete derivations of the NS equations from molecular models: however 
all attempts (which achieve the result under reasonable extra assumptions) deal with limiting 
regimes implying restrictions on initial data and involving a length scaling of 0(<5 _1 ), a time 
scaling of 0(5~ 2 ), (hence) a velocity scaling 0(5) and become exact in the limit as 5 — > (6 
being a dimensionless scaling parameter). 

Here we shall assume that the NS equations can be also obtained from a molecular model 
under a suitable scaling of space and time variables. We shall therefore consider microscopic 
initial data with velocity and position distributions with average fields (of density p(x), 
of kinetic energy, i.e. temperature, T(x) and velocity u(x)) consistent with initial values 
corresponding to a continuum. And we shall suppose that they evolve so that average 
velocity, density, kinetic energy satisfy NS with good approximation, and exactly in the 
limit in which some scaling parameter 5^0. The temperature field is identified with the 
average kinetic energy per particle, with velocities measured with respect to the average 
velocity u. 

Physically 5 is a parameter measuring "how far from a continuum the microscopic struc- 
ture is"; it can be identified with the ratio between the molecular free path and the length 
scale of the variation of the macroscopic velocity and temperature fields. 




(5) 



differ from this by -| J dCo ds,. 
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Therefore in the limit 5^0 each volume element will contain an infinite number of 
particles and fluctuations will be suppressed; however the average entropy creation will be 
defined and, by Eq.(3), be 



where k is the thermal conductivity, n is the outer normal to <9C and the average is intended 
over a time scale long compared to the microscopic time evolution but macro scopically short. 

The subscript ^ indicates that the time average (e)^ of Eq.(5) coincides with the average 
with respect to a probability distribution ix that is defined on the microscopic phase space and 
that we assume to be identified with a kind of evolving SRB distribution representing local 
stationarity. I.e. we assume the validity of a "Chaotic hypothesis", [1], not only globally 
and asymptotically in time but also locally and over times long compared to the microscopic 
time evolution but macroscopically short. This means that ^ may change with time over 
macroscopic time scales and is, therefore, in general different from the true SRB distribution 
that would describe the system over time scales very long with respect to macroscopic times. 

Hence suppression of fluctuations will not mean that the //-averages, over times long with 

def 

respect to microscopic scales, defining u(x),T(x) or the average U(x) = (U) (x),x E C , 
of the microscopic energy density U or the average energy per unit surface of the external 
thermostats U ext (£) = {U ext ) ($), £ € dC , over such scales will not continue to vary over 
macroscopic times. 

In other words we are saying that there are two basic time scales: at times long with 
respect to the microcopic time scale (determined by the molecular free path) the volume 
elements of the fluid are in a stationary state and the distribution fi used to compute 
the averages can be thought as being locally stationary on times short compared to the 
macroscopic times. As time evolves on the macroscopic time scale the distribution n also 
changes: to average, asymptotically, to the true SRB distribution for the global stationary 
state. 

Eq.(6) is the expression corresponding to Eq.(3) derived from molecular dynamics and it 
must be compared, for compatibility, with the familiar expressions for the entropy creation 
rate in systems described by macrosopic continua equations, [3]. 

We consider here a viscous and thermally conducting NS-fluid in local equilibrium. Thus 




(6) 
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the local equilibrium entropy density s depends on the local temperature and density: s = 
s(T,p). Let tJ be the stress tensor = (diUj + djui) in terms of the velocity field u, rj 
be the dynamical viscosity and U(x) = (U) (x) be the internal energy density (sum of the 
potential energy density and the kinetic energy density evaluated with respect to the average 
velocity u). Then the NS equation are, [4, p. 2, 3, 6, 18], 



(1) d tP + d-(pu) = 

(2) d t u + u-du = -±dp+2Au + g (7) 

(3) d t U + d- (uU) =r]r'du + kAT -pd-u 

(4) T(d t s + d- (us)) = nr'du +kAT 

here g is a (nonconservative) external force generating the fluid motion and p is the physical 
pressure. The conditions at the boundary of the fluid container C will be time independent, 
T = T(£) and u = (no slip boundary). 

As mentioned, Eq.(7) are macroscopic equations that can be valid only in some limiting 
regime. Given a system of particles with short range pair interactions let 8 be the mentioned 
natural dimensionless scaling parameter; then a typical conjecture is: for suitably restricted 
and close to local equilibrium initial data (see [4, p. 21] for examples) on time scales ofO(S~ 2 ) 
and space scales 0(S~ 1 ) the evolution of p, u, T, U, s follows the incompressible NS equation, 
[4, p.30]. 

The classical entropy creation rate in nonequilibrium thermodynamics of an incompress- 
ible fluid is, [4, p. 6], 



f dT 1 

k B £dassic = J c (« (^r) 2 + V Tfl'du) dx. (8) 

By integration by parts and use of the first and fourth of Eq.(7), k B e c i assic becomes, if 

def 

S = J c s dx is the total thermodynamic entropy of the fluid, 



J [-KdT ■ dT' 1 + r]^r' du)dx = 

(kAT + r)T'du) 
Jc T 
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n-dT , 
k — — — dst 

dC T 



dx = 



(9) 



id 
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k — — — ds£ + b + / u-osdx = 
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ie. this stiZZ /eads to the expression Eq.(6), "local on the boundary" or "localized at the 
contact between system and thermostats" , since u • n = by the boundary conditions, plus 
the time derivative of the total "thermodynamic entropy" of the fluid. 

Remarks: (i) An identical analysis can be performed for Rayleigh's convection model, widely 
used to test ideas on turbulence since [5]: the result is the same because the extra term that 
would appear in Eq.(9), see [4, p. 47], would be proportional to f Co u z dx which vanishes 
because the motion has no net momentum in the z direction. 

(ii) It should be noted that in the limit 8 — > 0, i.e. when the NS equations are expected to 
become rigorously exact, the Eq.(8) simplifies: only the first term in r.h.s. remains because 
the velocity u scales as 0(5), [4, p. 26]. 

3. Incompressible continua viewpoint 

The above analysis leads to a further natural question: whether the phase space contrac- 
tion and the entropy creation rate can be computed if we consider an incompressible fluid 
in local equilibrium and observed on a short but macroscopic time scale and imagine, as it 
is tempting to do and often done, [3, 4], its small macroscopic volume elements as a col- 
lection of small thermostats in contact with reservoirs consisting in the neighboring volume 
elements: if we apply the general theory of Sec.l would the results be consistent with the 
ones in Sec. 2? 

In establishing a comparison it should, however, be noted that the volume elements 
E = dx are not separated by walls, hence they can exchange particles, and they also move. 
In order to be able to treat volume elements as systems on their own we imagine that their 
size is A with A macroscopically small but microscopically large: certainly such length scale 
A is <C L ~ (f;fg) -1 where L is the macroscopic scale of the container C . 

Furthermore we have to assume that molecules diffuse, in a characteristic evolution time, 
over a distance <C A. The diffusion coefficient is D = 0{ ^£' ) with v the average speed, 
v = 0(v sound), m the mass of the molecules, r their radius and p the numerical density, and 
a characteristic time scale is 9 = f^. The distance traveled by diffusion in the latter time 
scale is (DO)* (~ KT 2 cm in air at normal conditions). 
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We assume, as above, that the local quantities, velocity field and temperature field, 
u(x),T(x), as well as the energy fields {U) ^(x) = U(x),x E C , and (U ext )^) = U ext (£),£ E 
dCo, evolve on a time scale much slower than the microscopic time scale 9 and can be 
considered constant on the latter time scale. The expression Eq.(8) should be regarded as 
an average over a long microscopic time (but over a short macroscopic time). 

Suppose that the conditions allow us to consider a volume element in a fluid as a ther- 
mostated system in a stationary state in contact with thermostats made of the neighboring 
elements. I.e. suppose that the quantity S introduced after Eq.(7) is small and the diffusion 
across the elements boundaries is not important enough to make the identity of the volume 
elements ill defined (i.e. 5 <C A <C L). Then we can apply the analysis leading from Eq.(l) 
to Eq.(3) and conclude that up to a total derivative R the phase space contraction of the 
total system, i.e. fluid plus thermostats, is (see also remark (ii) in Seel) 



<*) = ^rf- no) 

E E' B E' 

where Qe,e' is the amount of work that the particles in a given volume element E perform 
over the neighboring elements E', see Eq.(3). The quantity e(x) will be the entropy creation 
rate defined by regarding the volume elements as small thermostats in stationary state. 

If the average heat current is —K,dT and the element E is imagined with the bases 
orthogonal to the gradient of T, the average contribution to Qe,e' for E' adjacent to the 
upper base of E is —KdT-nX 2 (A 2 = area of the base) and it is opposite to the contribution 
from the lower base; therefore the quantity J2e' r E f nas avera S e 



-k&T • n\ 2 (-^- — -^-) = —ndT • dT~ l \\ (11) 

if T± are the temperatures at the two bases; hence summing over E: fcs(e) = J Co d(^) ■ 
(-K&T) dx which can be written in the more common and familiar form k J Cq (^-) 2 dx. Then 
if d t T = 0, AT = and u = Eq.(8), therefore Eq.(6), follows by partial integration. 

More generally in presence of time dependence and non vanishing velocity field there will 
be an extra amount of energy transfered to elements adjacent to E and due to diffusion 
across the bases: it can be evaluated in the same way as above to be r/r' ■ n (u + — u_) A 2 if 
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r ' is the stress tensor and it changes k B {s) to Eq.(8), hence to Eq.(6) plus S (see Eq.(9)) 
if account is taken of the fourth of Eq.(7). 

The above viewpoint therefore also leads to the classical expression for the entropy cre- 
ation rate. 

4. Phase space contraction 

Finally it is interesting to remark that not only the entropy creation rate but also the 
phase space contraction can be computed along the above lines for macroscopic continua. 
Regarding each volume element E as a thermostated system in a stationary state with a 
fixed temperature, the average phase space contraction is deduced from Eq.(3) taking into 
account the R contribution. 

Hence, setting ks e(x) d = (^f ) 2 + rj^r'du^J, the average phase space contraction is 

Jco ( e fe) — ksTix) ) ^— ~~ fdc fcflT(g) ^ s £> wnere U(x) = (E/) M (x) denotes the average potential 
energy density while with U ext (t) = (U ext ) (£) we denote the average potential energy density 
per unit surface of the thermostats. 

Remark that the time derivative of the potential energy density U (x) could be replaced 
by the time derivative of the total energy density at x: in fact the kinetic energies do not 
contribute to the total derivative of the energy density since they have been supposed to 
be constant because each volume element is regarded to have a constant kinetic energy, i. e. 
a well defined temperature, and in the above relations the energies appear through their 
time derivative). Hence energy conservation, U(x) = rjr' • du + kAT, see (3) in Eq.(7), for 
incompressible fluids and partial integration of the contribution k(=^) 2 to e(x), see Eq.(8), 
in the integral J Cq (e(x) — k^r{ x ) ) d% leave us with a boundary term which is just Eq.(6) 
minus J dC(j j^y^j ds^ as it could be guessed from the expression for a preceding Eq.(3). 

Therefore even if we regard the fluid volume elements as thermostated systems in sta- 
tionary nonequilibrium we are led to the expected relation between average entropy creation 
(e) and average phase space divergence (<r) , namely (if x denotes the fields determining 
the state of the fluid) 



<"><•=<«>*-/.* Mr- (12) 

Analysis of compressible fluids, unfortunately more difficult, should also be attempted: 
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the first difficulty will be, of course, that it is not clear under which scaling the compressible 
NS equations should hold as a reasonable approximation. 

5. Conclusion 

The conclusion is that the entropy creation rate evaluated from the microscopic model, 
Eq.(l), and the microscopic definition of heat ceded to a thermostat (see Eq.(2)) 

(1) differs "by a total derivative", R, from the microscopic phase space contraction (see 
Eq.(3)), 

(2) if the system can be regarded as a NS continuum in stationary local equilibrium, its 
average coincides with the classical entropy creation rate up to a total derivative, S, see 
Eq.(9), 

(3) in the latter situation its average coincides with the entropy creation obtained by re- 
garding the continuum as constituted by volume elements each of which is a thermostat in 
a stationary state exchanging heat with the neighboring volume elements, see Eq.(8), and 
differs by a total derivative from the average phase space contraction evaluated by regarding 
the continuum in the same way, see Eq.(12). 
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